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GRAVITY MODEL IMPROVEMENT USING 
GEOS- 3 (GEM 9 & 10) 


ABSTRACT 

The spocebome altimeter missions of GEOS- 3 (50 cm occur acy) and the 
future SEASAT (10 cm occur ocy) require precise knowledge of the radial position 
of the spacecraft to be most effective. Though errors in previous gravity models 
have produced large uncertainties in the orbital position of GEOS- 3, significant 
improvement has been obtained with new geopotential solutions, Goddard Earth 
Models 9 and 10. Using least squares collocation GEM 9 was derived by 
combining laser data from GEOS-3, LAGEOS and Storiette, S-Band measure- 
ments on LANDSAT I, together with data from 26 other satellites used in 
previous solutions. GEM 10 is a combination solution containing a global set of 
surface gravity anomalies along with the data in GEM 9. Radial errors of GEOS- 
3 for 5 day arcs have been reduced from about 5 m to I m based upon orbital 
intercomparisons, station navigations and analyses employing crossover points 
from passes of altimetry. 

The use of collocation has permitted GEM 9 to be a larger field than 
previous derived satellite models, GEM 9 having harmonics complete to 20 x 20 
with selected higher degree terms. The satellite data set has approximately 
840,000 observations, of which 200,000 are laser ranges taken on 9 satellites 
equipped with retroreflectors. GEM 10 is complete to 22 x 22 with selected 
higher degree terms out to degree and order 30 amounting to a total of 592 
coefficients. Comparisons with surface gravity and altimeter data indicate a 
substantial improvement in GEM 9 over previous satellite solutions; GEM 9 is in 
even closer agreement with surface data than the previously published GEM 6 
solution which contained surface gravity. In particular the free air gravity 
anomalies calculated from GEM 9 and a surface gravity solution by Rapp (1977) 
are in excellent agreement for the high degree terms (13 < l< 22). 
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The mass constant of the Earth, GM, has been estimated from the laser 
3 2 

data as 398600.64 + .02 km /sec , a value which is principally determined from 
LAGEOS. The speed of light used was 299792.5 km/sec. Geocentric station 
positions were determined for approximately 150 stations in GEM 10. These 
station coordinates, their mean sea level heights and altimetry data provide an 
estimate for the mean radius of the earth of a fi = 6378140 ♦ I m. Accuracy 
estimates derived for the potential coefficients have been verified with 
independent data sets. These produce commission errors in geoid heights of 1.9 m 
and 1.5 m (global RMS values) respectively for GEM 9 and 10. 
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SECTION I 


INTRODUCTION 



I. INTRODUCTION 


The Earth and Ocean Dynamics Applications Program (EODAP) of the 
National Aeronautics and Space Administration calls for knowledge of the global 
geoid to sub-meter levels of accuracy. While final realization of these goals will 
rely strongly on GEOS-3 and SEASAT altimetry, progress continues to be made 
toward comprehensive gravity modeling using conventional satellite tracking 
systems and surface gravimetry. 

At Goddard Space Flight Center (GSFC) the emphasis has been on using 
as much of the precise satellite data as possible. Precise laser tracking such as 
the International Satellite Geodesy Experiment (ISAGEX, Brachet, 1970) laser 
systems of 1970 have yielded a substantial improvement in geopotential 
sensitivity and accuracy over the last few years (e.g., Wagner, et al 1977). 
However, the accuracy of the GSFC, SAO and French laser systems (with 5 cm 
noise levels for GSFC systems now deployed) on the new GEOS-3 and Storiette 
orbits could not be realized without continued geopotential improvement. In the 
case of GEOS-3, effective use of the altimeter data required a very significant 
improvement in radial orbit determination accuracies beyond the capabilities of 
existing gravity models. Improvement of GEOS-3 orbit determination by 
reduction of geopotential uncertainties was a major objective of Goddard Earth 
Models (GEM) 9 and 10. 

GEM 9 is a gravity model based solely on optical, laser, and electronic 
observations taken on 31 satellites. GEM 10 combines the GEM 9 satellite data 
with surface gravimetry. (GEM 10 and other solutions which are derived from 
both satellite and surface observations (e.g., SAO 4.3, GRIM 2) are referred to as 
"combination" solutions. ) 



GE.Vi 9 and 10 incorporate a number of significant changes in technique 
over previous GEM solutions. The extension of the GEM 9 satellite solution to 20 
x 20 (complete in degree and order) was accomplished through the use of least 
squares collocation (Moritz, 1972). This technioue is discussed in Section 3.2 and 
is also used by King-Hele (1974) and Anderle (private communication, 1977) in 
their gravity work. The adjustment of the earth's mass (GM) is another 
advancement (Section 3.3). A significant improvement was obtained in GEM 10 
by now including the truncation of the gravity field (as well as the accuracy of 
the data) as an error source in weighting the gravimetry observations (Section 
3.4). GEM 9 and 10 will be used as the base fields for other solutions being 
planned which will extensively use the oltimeter data available on GEOS-3. 

Many of the data systems for the GEOS-3 mission have been used to 
evaluate the GEM 9 and 10 models. While the satellite-to-satellite doppler relay 
(SSt) ar»d the altimeter ranging data have not been included in these latest GEM 
models, these data have been used to assess the overall global improvement of 
the models. The laser tracking has aLo been used to test the models. These 
studies are included within this report end provide a strong demonstration of the 
high level of accuracy which has been achieved in GEM 9 and 10. 
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SECTION 2. 


DATA 




2. DATA 


2.1 SATELLITE TRACKING DATA 

A brief summary of the 840,000 satellite tracking measurements utilized 
in GEM 9, is .liven in Table I. The main feature of the data in the new solution is 
the large amount of laser data employed totaling about 200,000 observations. 
Because of the sensitivity of the laser system to satellite perturbations (down to 
5 cm), contributions of the laser observations have been computed complete 
through degree and order 22 for the harmonics, whereas the harmonics were 
computed complete only through degree and order 16 for the other types of data. 
The ISAGEX laser data have been used in previous solutions, but in these, the 
harmonics were computed complete only through degree 16. 

A description o f lie satellites employed and their data distribution is 
given in Tables 2, 3A, and 3B. These tables respectively describe (2) satellite 
orbital characteristics and types of data employed, (3a) the distribution of data 
on satellite arcs containing optical data only, and (3b) the distribution of data on 
satellite arcs containing a variety of tracking systems consisting of electronic, 
laser, and additional optical observations. 

Characteristics of the data among the various tracking systems are 
summarized in Table 4. Summaries by tracking network consist of the number of 
stations, observations, and satellites observed including accuracies and weights 
used for sigmas of the data in the solution. There are 561,000 measurements 
which have been used previously in GEM 7 and these are distributed among 9 
different tracking networks. The table also shows the data which are unique to 
GEM 9, totaling 278,400 observations for Laser, S-Band, and NWL Doppler 
tracking systems. 
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TABLE 1. 

GEM 9 SATELLITE TRACKING DATA 


NO. HARMONICS 


TYPE 

NO. OBS. 

SATELLITES 

NO. ARCS 

(Compfeta) 

OPTICAL 

150.000 

24 

287 

16 x 16 

ELECTRONIC 

477.000 

11 

97 

16 x 16 

LASER 

213,000 

9 

127 

22x22 


LASER DATA DISTRIBUTION 


GEOS— 3 

94,000 obs. 

38 arcs 

STARLETTE 

28,000 

26 

LAGEOS* 

25,000 

11 

CD 

m 

1 

o 

3,000 

4 

7 ISAGEX SATELLITES 
(BE-B. BE-C. D1-C. 
D1— D, GEOS— 1, 

GEOS— 2, PEOLE) 

63,000 

48 


•LAGEOS USED FOR ESTIMATING GM AND STATIONS ONLY (SEE SECTION 3.3). 
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TABLE 2. SATELLITE ORBITAL CHARACTERISTICS USED IN GEM 9 AND 10 


PRIMARY 

MEAN RESONANT 

SATELLITE A (KILO- I MOTION PERIOD 


NAME 

METERS) 

E 

(DEGREES) 

(REV/DAY) 

DAYS 

DATA TYPE 

AGENA-RB* 

7297. 

0.0010 

69.91 

13.92 

5.0 

o 

ANNA— IB 

7501. 

0.0082 

50.12 

13.37 

4.8 

O. RR 

BE-B 

7354. 

0.0135 

79.69 

13.76 

3.0 

L, RR, O 

BE-C 

7507. 

0.0257 

41.19 

13.35 

5.6 

L, RR, O 

COURIER 

7469. 

0.0161 

28 31 

13.46 

3.8 

O 

Dl-C 

7341. 

0.0532 

39.97 

13.81 

2.5 

L.O 

Dl-D 

7622. 

0.0848 

39.46 

13.05 

8.4 

L, O 

ECHO— 1RB 

7966. 

0.0118 

47.21 

12.21 

11.9 

O 

GEOS— 1 

8075. 

0.0719 

59.39 

11.96 

7.0 

L, RR, O 

GEOS— 2 

7711. 

0 0330 

105.79 

12.82 

5.7 

L, R, RR, O 

GEOS— 3 

7226. 

0.0008 

114.98 

14.13 

4.5 

L 

GRS 

7239. 

0.0598 

49.76 

14.10 

10.7 

O 

INJUN 

7316. 

0.0079 

66.82 

13.87 

3.8 

O 

LANDSAT-1 

7286. 

0.0013 

99.10 

13.99 

18.0 

RR 

LAGEOS 

12273. 

0.0038 

109.85 

6.39 

2.7 

L 

MIDAS-4 

9995. 

0.0112 

95.83 

8.69 

3.0 

O 

OGO-2 

7341. 

0.0752 

87.37 

13.79 

3.8 

O 

OSCAR -7 

7411. 

0.0224 

89.70 

13.60 

2.2 

o 

OVI-2 

8317. 

0.0184 

144.27 

11.45 

2.2 

0 

PEOLE 

7006. 

0.0164 

15.01 

14.82 

2.1 

L, M 

SAS 

6923. 

0.0035 

3.04 

15.09 

4.6 

M 

SECOR-5 

8151. 

0.0793 

69.22 

11.79 

3.4 

O 

STARLETTE 

7331. 

0.0204 

49.80 

13.83 

2.8 

L 

TELSTAR 

9669. 

0.2429 

44.79 

9.13 

14.9 

O 

TIROS— 9 

8024. 

0.1173 

96.41 

12.07 

19.5 

M 

TRANSIT-4A 

7322. 

0.0076 

66.82 

13.85 

3.5 

O 

VANGUARD-2RB 

8496. 

0.1832 

32.92 

11.09 

294.3 

O 

VANGUARD-2 

8298. 

0.1641 

32.89 

11.49 

2.7 

O 

VANGUARD-3 

8508. 

0.1901 

33.34 

11.07 

187.6 

O 

5BN-2 

7462. 

0.0058 

89.95 

13.46 

2.4 

O 


*RB ^ Rocket Body 

## L - Laser Range, R — Range, RR — Range Rate, O - Optical, M — Mini track 
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TABLE 3A. DISTRIBUTION OF DATA FOR SATELLITE ARCS 
USING OPTICAL DATA ONLY 

287 WEEKLY OPT. ARCS (PRIMARILY SAO BAKER-NUNN) 


SATELLITE 

NAME 

& SATELLITE 

ID 

AGENA-RB 

640011 

ANNA— IB 

620601 

BE-B 

640841 

BE-C 

650321 

COURIER— IB 

600131 

Dl-C 

670111 

O 

1 

Q 

670141 

ECHO— IRB 

600092 

GEOS— 1 

650891 

GEOS— II 

680021 

GRS 

630261 

INJUN-1 

610162 

MIDAS-4 

610281 

OGO-2 

650811 

OSCAR-7 

660051 

OVI-2 

650781 

SECOR-5 

650631 

TELSTAR-1 

620291 

TRANSIT-4A 

610151 

VANGUARD-2RB 

590012 

VANGUARD-2 

590011 

VANGUARD-3 

590071 

5BN-2 

630492 


NO. NO. 

ARCS OBS. 


7 1005 

40 4183 

4 469 

22 4947 

12 3375 

4 902 

9 6386 

18 2240 

28 40855* 

24 25315* 

5 369 

9 768 

20 14879 

7 461 

4 1780 

4 910 

4 290 

16 1946 

14 1316 

11 379 

5 615 

15 996 

5 355 

TOTALS 287 114700 


•MOTS/SPEOPTSOBS.: GEOS-1 - 22100, GEOS-II - 22000 PLUS 2100 OBS. FROM 
INTERNATIONAL CAMERAS. 
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TABLE 3B DISTRIBUTION OF DATA FOR SATELLITE ARCS USING ELECTRONIC AND 
LASER TRACKING SYSTEMS OR COMBINED SYSTEMS INCLUDING CAMERAS 
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DAY LONG ORBITAL ARCS 
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2.2 SURFACE GRAVITY DATA 


A set of 1654 equal area 5° mean gravity anomalies (Rapp, 1977), have 
been used along with the satellite tracking data in ^He combination solution GEM 
10. The data is based upon approximately 38,000 1° mean gravity anomalies 
(Figure 13). Accuracy estimates for the 5° mean anomalies are depicted in 
Figure I. Of the 1654 5° mean anomalies, 1507 were based directly on the 1° 
anorriHies whrft* the remaning 147 5° means were obtained by interpolation. The 
distribution of the number (N) of 1° anomalies within a 5° block is shown in Table 
5 along with accuracy estimates of the 5° means. Only 625 of the 5° blocks 
contain a full set of 1° mean (observed) anomalies. 


II 



Figure t, Rapp 5° {equal area) Anomalies and Their Uncertainties 
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TABLE 5. AVERAGE ACCURACY OF 5° MEAN ANOMALIES COMPARED TO 
THE NUMBER (N) OF 1° ANOMALIES WITHIN THE 5° BLOCK 


N 

NUMBER OF 
5° MEANS 

AVERAGE 

ACCURACY 

(MGALS) 

25 

625 

2.5 

20-24 

310 

3.5 

15 - 19 

177 

5.3 

10 - 14 

151 

7.2 

5-9 

144 

10.0 

1 -4 

100 

14.0 

0 

147* 

17.0 

TOTAL 

1654 



•INTERPOLATED FROM NEIGHBORING 5° ANOMALIES. 


13 




SECTION 3. 


MODELING TECHNIQUES AND RESULTS 
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3. MODELING TECHNIQUES AND RESULTS 


The basic modeling techniques employed for the orbital, geopotential, 
and station solutions are given in detail by Lerch et al, 1974. In this section we 
present the GEM 9 and 10 solutions for the potential coefficients and station 
coordinates along with a discussion of the new techniques employed. We also 
present and discuss solutions for fundamental geodetic reference parameters: 
the mean radius of the ear+h (a ), the gravitational constant (GM), and mean 
equatorial gravity (g e ). 


3.1 GEOPOTENTIAL 


The gravitational potential was modeled in terms of spherical harmonics 
as follows: 

•3 

l 


- !'• EZC J ) 

r / t * „ A / 


+ S^f sin mX 


/ -2 m=0 

I 


% <sin 0) 


C a cos m\ 
£ m 


: m 


i 


where GM is the earth's gravitational constant including the atmosphere, a^ is 
the earth's mean equatorial radius, is the fully normalized associated 

Legendre function of degree t and order m (e.g., Kaula, 1966, p. 7) and r,o, X 
are the distance from the center of mass, latitude and longitude. The normalized 
potential coefficients (C^ , tor GEM 9 and 10 are listed in Tables 6 and 7. 
Using these potential models in Brun's formula (Heiskanen and Moritz, 1967) 
geoids are computed and presented in Figures 2 and 3. 
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GEOPOTENTIAL COEFFICIENTS 



TABLE 6. GODDARD EARTH MODEL 9 
NORMALIZED COEFFICIENTS (UNITS OF 10“' 
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0.07315 0.05121 9 2 0.02445 -0.02467 10 2 -0.08429 -0.00376 

0.03947 -0.09951 12 2 0.00390 -0.00931 13 2 0.02214 -C. 03804 

-0.03632 0.04590 15 2 0.00864 -0.02976 16 2 -0.00624 0.01666 

-0.00985 0.C1244 18 2 0.00649 0.02368 19 2 0.00877 -0.00714 

-0.00395 0 00673 21 7 0.01467 0.00653 22 2 -0.00147 -0.00964 



TABLE 6. GODDARD EARTH MODEL 9 (continued) 
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TARLE 6. GODDARD EARTH MODEL 9 (continued) 
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-0.00851 -0.00934 19 14 -0.00552 -0.01273 20 14 0.01232 -0.01019 

0.01955 0.01031 22 14 0.00954 0.00665 23 14 0.00881 -0.00568 

-0.01659 0.00317 25 14 -0.02360 0.01669 26 14 0.00688 0.00148 



TABLE 6 . GODDARD EARTH MODEL (continued) 
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TABLE 7 . GODDARD EARTH MODEL 10 
NORMALIZED COEFFICIENTS (UNITS OF 10 -1 
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TABLE 7. GODDARD EARTH MODEL 10 (continued) 
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3.2 


COLLOCATION 


The major innovation in GEM 9 over previous Goddard Earth Models was 
the use of "least squares collocation" (Moritz, 1972), which allowed the extension 
of the satellite field to 20 x 20. In this procedure we employed an approach 
similar to that of Rapp (1973; eq. 13). Conventional least squares simply 
minimizes the observation residuals (noise); the results of such an approach are 
described graphically in Figure 4 for gravity model recovery. In this figure the 
solution without any constraints (simple least squares) diverges at high degrees 
from the independent surface gravity data used to test it. The high correlation 
between certain high degree and order coefficients is the problem which causes 
an excessive adjustment of the coefficients in the solution. Least squares 
collocation essentially minimizes both the signal (e.g., harmonic coefficients) and 
the noise (observation residuals), thus controlling the excessive adjustment. 
First, we present the technique. The result of its application for GEM 9 and the 
result of other tests shown in Figure 4 are then discussed. 

The principle of collocation is to minimize 

r T Wr + s T Ws £ Q (I) 

with respect to the unknowns x, where 

x - geopotential, station and orbit parameters 

r - satellite observation residuals 

W - diagonal weight matrix for satellite observation residuals 

s - signal, harmonic (potential) coefficients representing a 
subset of x 

2 

W - diagonal weight matrix with elements 1/^ 
where T s(/,m) = ICf'V/^ 
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MEAN SQUARE RESIDUAL: MGAL 


FIGURE 4. COMPARISON OF TRUNCATED GEM 9 SATELLITE DERIVED FIELDS 
USING DIFFERING LEVELS OF COLLOCATION WITH SURFACE GRAVITY DATA 
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Let s represent a subset of the potential coefficients with the partition 



y 


V 

X = 


- 



s 


-IS 


( 2 ) 


and using the linearized forms from Taylor's series 


r = r -A.Ay-B.As 
o ' 


s = s„ + -As, 
o 


(where A and B are matrices 
of partial derivatives) 


then minimizing Q in (I) above gives the normal equations 


A T WA 

a t wb 




A T Wr 

o 

b t wa 

b t wb + w 


-AS 


B T Wr -Ws 
L 0 O J 


(3) 



Allowing for a scale tactor w to adjust the relative weighting between W and W 
above, we have 


W = w W (W is the formal weight matrix) 
o o 3 


; f 2/,2 

w = f /f 


(5) 

( 6 ) 


" = f,7 o 



where f is an estimate for scaling up the standard deviaiions (^ , of the potential 
coefficients implicit in the satellite normal equations and f is a corresponding 
estimate for scaling the rms size coefficient ( * ) os given by Kaula's rule. Based 
upon the size of the coefficients and the scaling of their standard errors ir, GFM 
7, we used f -yf \ /2 and f = giving w = .05 in GLM ?. 
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TABLE 8. RATIO OF DIAGONAL TERMS (d> OF THE SATELLITE NORMAL MATRIX 
IN GEM 9 TO THE DIAGONAL TERMS (d) OF THE SIGNAL MATRIX 


d 


do -5 // 2 ) 


— 4 . 


FOR DEGREE 


l 


m* 

1 = 12 

1 = 16 

1 = 20 

0 

1,000,000 

250,000 

20,000 

1 

80,000 

6,000 

630 

2 

13,000 

2,500 

630 

3 

10,000 

2,000 

200 

4 

5,000 

1,000 

200 

5 

8,000 

800 

160 

6 

2,500 

1,300 

160 

7 

10,000 

300 

160 

8 

4,000 

1,300 

80 

9 

13,000 

800 

250 

10 

20,000 

1,600 

80 

11 

20,000 

2,500 

40C 

12 

100,000 

25,000 

4,000 

13 


40,000 

16,000 

14 


630,000 

63,000 

15 


100,000 

2,500 

16 


600 

1,600 

17 



800 

18 



310 

19 



80 

20 



25 


# C and S tesseral terms are essentially the same. 
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Since the signal matrix contains only diagonal terms which were added to 
the data matrix it is interesting to compaie their relative sizes* As seen from 
Table 8, the satellite unsealed normal equations (B^W q B) have considerably 
larger diagonal terms (even out to degree 20) than the signal matrix W, which is 
still true even after w = .05 is applied. This demonstrates that our application of 
collocation can have a significant effect only by indirectly controlling ill- 
conditioned vectors (correlation effects) in the system. 

Collocation was applied to the coefficient subset for degree l > 12 
except for resonant terms of order 12, 13, end 14. 

Figure k shows the improvement when collocation is applied to terms 

above degree 12 compared to those when applied above degree 16. The former 

solution was chosen as GEM 9. Interestingly, when collocation was applied to 

terms above degree 8, the results were almost the same indicating that this 

method was unnecessary for the lower part of the recovered geopotential. The 

results for GEM 7 in Figure 4 show that simple least squares can provide a 

reasonable satellite solution complete to 16 x 16. A solution similar to GEM 7 

(16 x 16, no collocation) was obtained using the GEM 9 data, giving results 

2 

comparable to GEM 7 within one mgal. Hence, without controlling matrix ill- 
conditioning it is unlikely that worldwide geopotential improvement would have 
resulted using the new data. GEM !0 was gIso derived using the collocation 
technique applied to the coefficients above degree 12. GEM 10 is complete to 22 
x 22. H is important to note (as mentioned above) that collocation was not 
app^ed to the resonance terms (m = 12, 13 and 14). 


3.3 DETERMINATION OF GM 

The simultaneous determination of GM with the geopotential and station 
positions wes performed. Table 9 describes additional tests which were made to 
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TABLE 9. 

GM DERIVED FROM SATELLITE LASER DATA 


SATELLITE 

GM (km^/sec^) 

NO. OF 
5-DAY ARCS 

NO. OF OBS. 

LAGEOS 

398600.64 

11 

25,000 

STARLETTE 

0.70 

26 

28.000 

GEOS— 3 

0.84 

38 

94,000 

COMBINED 

398600.64 




SUBSETS 

GM 

NO. OF ARCS 

NO. OF OBS. 

STARLETTE 

398600.44 

9 

12,000 

STARLETTE 

0.87 

9 

10,000 

STARLETTE 

0.73 

8 

6,000 

GEOS— 3 

0.85 

18 

33,000 

GEOS— 3 

0.92 

16 

54,000 

LAGEOS 

0.64 

5 

16.000 

LAGEOS 

0.65 

6 

9 000 


LAGEOS SUBSETS 

DETERMINED 

FORMAL 

STANDARD 

ERROR 

NO. OF 


GM (km 3 /sec 2 ) 

(km 3 /sec^) 

OBS. 

LAGEOS INDIVIDUAL ELEVEN ARCS 

1 

.635 

.014 

2.037 

2 

.567 

.054 

651 

3 

.538 

.020 

2,830 

4 

.682 

.095 

1,167 

5 

.664 

.011 

2.203 

6 

.687 

.020 

2.037 

7* 

.493 

.431 

1,676 

8 

.641 

.032 

1,492 

9 

.602 

.016 

6.021 

10 

.647 

.019 

1.481 

11 

.629 

.020 

3.634 

ARCS 1 THROUGH 6 <6 arcs) 

.648 

.007 

10.880 

ARCS 7 THROUGH 11 <5 arcs) 

.647 

009 

14,259 

WORST CASE (from LAGEOS) 

6 HIGH GM ARCS 

.661 

.007 

11,960 

5 LOW GM ARCS 

.611 

.008 

13,179 

ALL DATA 

398600.638 

.005 

25,139 


# Inadvertantly, this was only a 10 hour arc. 
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evaluate solutions for GM where only station coordinates and orbit parameters 
are solved simultaneously with GM. The speed of light used was 299792.5 
km/sec. The estimation of GM was exclusively from laser tracking data. The 
presence of LAGEOS dominated the determination, and the GM results obtained 
were repeated using the LAGEOS data by itself. The value of GM is 398600.64 
km /sec for all 1 1 LAGEOS arcs. LAGEOS not only dominated the combination 
solution, but in the subset solutions it also gave much more consistent results 
than either GEOS-3 or Storiette. This is because the high altitude of LAGEOS 
provides good geometry and dynamics for estimating GM with separability for 
station coordinates. The individual LAGEOS arcs shown in Table 9 were 
recombined taking the highest 6 determined values for one solution and the 
remaining 5 lowest values in a second solution. These two "worst case" solutions 
were both within .03 of the above value of GM, whereas a typical set of arcs 
(first 6 and last 5) are within .003 of this value. Based upon these results and the 
formal uncertainty of .005 for the total solution value (398600.64), the value of 
.02 was selected as a conservative estimate of the uncertainty for GM. 


3.4 MODIFIED TREATMENT OF THE SURFACE GRAVITY DATA 
FOR INCLUSION INTO GEM 10 

Another major innovation over previous GEM solutions is in the 
treatment of the surface gravimetry. In GEM 10 the surface gravity data has 
less overall weight than in previous GEM combination solutions. For the GEM 10 
solution, an additional 5 mgal was added to each individual observation 
uncertainty. This 5 mgal uncertainty was used to represent the unmodeled 
truncation error for 5° mean anomalies when solving for a 22 x 22 field. This 
weighting scheme had the benefit of making the data quality more uniform over 
the globe than in previous models. The result was a solution which agreed with 
the gravimetry over the oceans about as well as over land. More importantly, 
agreement with worldwide altimetry was superior with this more uniform 
weighting (see later Section 4.4). 
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3.5 STATION COORDINATES AND GEODETIC REFERENCE PARAMETERS 


GEM 9 and 10 simultaneously determined the center of mass positions for 
146 tracking stations. These station coordinates for GEM 10 are presorted in 
Table 10. Table I I compares the GEM 10 station positions with those estimated 
by Marsh (1977) for the Calibration Area lasers. The geocentric station 
coordinate differences are seen to be about I m in these results. These results 
indicate the highly accurate laser station coordinates have been obtained, with 
uncertainties being significantly less than 3 m given in Lerch et al, 1974. 

Three methods were used to derive a mean value of the semi-major axis, 
a e > °f fbe earth's reference ellipsoid, all of which agree to within one meter of 
<* e = 6378140 m. These results made use of reference parameters such as GM, 
equatorial gravity (g^) and ellipsoidal flattening. These parameters are all 
compared in Section 3.5.4 with the set adopted by a special study group of the 
I AG in 1975. 


3.5.1 a Derived from the GEM 10 Station Coordinates and Their Mean 
e 

Sea Level Surveyed Heights 

Station coordinates and mean sea level heights from survey were used 
and gave - 6378139.9 + 1.5 m (Table 12). Subset solutions presented in Table 
12 for the different tracking systems all agree to within one meter of the mean 
value except for the Baker-Nunn sites which difler by 2.9 m. These results are 
based upon the following formula: 


a - a (reference) = 
e e 


£ 

r t 


h. - MSL H. - N. 

i i i 


L 
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TABLE 10. STATION COORDINATES OF GEM 10 (a. - 6378145m.. 1/1 - 296.2551 


STATION 


LA TI TOTE 


LONG 

I TUDE 

HEIGHT 

NAME 

NUMHEK 

m> 

MM 

ss.sss 

UGD 

MM 

SS.SSS 

METERS 

IBPOIN 

1021 

vs 

?5 

49.815 

282 

54 

49. 087 

-38.8 

1FTMYK 

1022 

?6 

32 

53.238 

278 

8 

4.616 

-25.6 

lOQMfcR 

1024 

-H 

? 3 

24.997 

l 36 

52 

15.624 

131.3 

1 GAT AC 

1026 

- * \ 

a 

58.433 

289 

19 

53. 929 

719.1 

1M0JAV 

1030 

*S 

19 

47 .928 

243 

5 

59.673 

893.9 

1 JORUR 

103 1 

-2S 

S3 

0 .858 

27 

42 

26.638 

1546.6 

INEWFL 

1032 

47 

44 

29.369 

307 

16 

46. 791 

70.1 

icntec 

1033 

64 

52 

18.279 

212 

9 

37. 760 

174.2 

IGFURK 

1034 

4 H 

1 

21 .325 

26? 

59 

20. 154 

224.2 

I VlNKEL 

1035 

si 

26 

45.911 

359 

18 

9.92 3 

ill. 5 

1ULASK 

1036 

S4 

68 

37.120 

212 

28 

31.481 

301.7 

1 rosmn 

1037 


12 

7.220 

2 77 

7 

41.811 

873.4 

1CRORL 

101ft 

* \*> 

3 7 

32.L20 

148 

57 

14. 996 

944.5 

iROSMA 

1042 

* r > 

1? 

7.1R5 

277 

7 

41.570 

872.3 

l T ANAN 

104 3 

-1 ) 

C 

31 .966 

47 

17 

59.667 

1371.6 

M AUG AH 

1122 

-11 

1 

14 .815 

47 

18 

11.514 

1382.3 

MAUf.AS 

1123 

-10 

1 

14 .306 

4 7 

18 

11.066 

1383.4 

RUSR AN 

1126 

35 

il 

45.502 

277 

7 

26.832 

837 .6 

ULASKK 

1 128 

64 


19.343 

212 

29 

13.468 

349.1 

CAR VON 

1152 

— ?4 

64 

10.5 70 

113 

42 

59.670 

12.0 

8CA3 

1 30 2 

\? 

21 

4.781 

295 

20 

31.728 

-16.9 

C Y 1 3 

l 304 

n 

45 

51 .725 

344 

/. 1 

58. 376 

193 .5 

CRP3 

l 30 8 

-?<* 

64 

23.279 

1 1 3 

4 3 

32.077 

16.5 

GURTkS 

1112 

22 

/ 

34.712 

200 

20 

5.478 

1158.1 

GUSA 

1314 

36 

20 

29.8 38 

24 3 

7 

37. 302 

924.4 

TEX3 

1316 

27 

>9 

17.723 

262 

37 

30. 991 

-34.3 

MAOH 

1323 

40 

77 

19.824 

355 

40 

53. 930 

826.8 

5 KM 3 

1324 

i 3 

1 8 

38 .148 

144 

44 

12.928 

136.6 

H SK 8 

l 325 

- 3S 

34 

59 .862 

148 

5H 

40.499 

113«.5 

GOSR 

1328 

35 

2n 

29 .879 

24 3 

7 

34. 870 

931.6 

MIL 3 

137 1 

?R 

30 

29.588 

279 

1R 

23.668 

-30.2 

AC93 

1375 

- 7 

5 7 

17.311 

345 

4C 

22.557 

551.5 

ETC3 

1377 

3* 

64 

54 .935 

283 

9 

2 0.064 

10. 0 

fcTCA 

1 19 1 

V> 

6 9 

54 .2 3 3 

26 3 

9 

29. 094 

16.3 

HUWAKI; 

200 1 

39 

4 

48 .4 78 

283 

6 

7.823 

119.3 

NEKMtX 

200 1 

3? 

16 

53.574 

25 3 

14 

53.601 

1165.5 

SANHtS 

2008 

-2 3 

1 3 

3.712 

314 

7 

4 >. 686 

595.6 

MISAWA 

2013 

4t) 

4 3 

14 .1 ?0 

1 4 1 

19 

51 . 729 

45.5 


37 



TABLE 10. (con'inuad) 


ANCHOR 

2014 

61 

17 

0 .052 

210 

10 

29. 747 

66.0 

T AEIJNA 

2017 

-14 

n 

50 . 1 o i 

189 

17 

3.411 

35.7 

TMULEG 

2018 

76 

32 

19 .93? 

291 

13 

53.664 

57.9 

MCMURU 

2019 

- 77 

50 

51 .687 

166 

40 

25.699 

-19.8 

AUSTIN 

2092 

30 

17 

13.432 

262 

16 

5.217 

156.9 

WAHIWA 

2100 

21 

31 

15.383 

202 

0 

10.710 

403.4 

L ACRES 

2103 

32 

16 

44 .153 

253 

14 

46.220 

1166.1 

LASHAM 

2106 

51 

11 

9.141 

353 

58 

25.656 

222.3 

APLKND 

2111 

39 

9 

48.273 

283 

6 

12. 323 

100.9 

SMITHL 

2112 

-34 

40 

26.262 

1 38 

39 

17. 124 

27.7 

PRETOR 

2115 

-25 

56 

48.1 70 

2P 

20 

52.011 

1597.0 

ASAMUA 

2117 

-14 

19 

50.257 

139 

17 

3. 354 

40.3 

S ANM I G 

2121 

14 

59 

16.402 

120 

4 

21.378 

58.7 

W AL UUP 

2203 

37 

51 

51 .793 

284 

29 

33.055 

-32.1 

CANTON 

2706 

-2 

47 

35.334 

188 

20 

4 . 890 

27.3 

MAHE 

2717 

-4 

*0 

13 .748 

55 

28 

46.830 

548.6 

A SC E MS 

2722 

-7 

58 

10 .006 

345 

35 

40.876 

92.1 

COCOS 

2723 

-12 

11 

44 .9 32 

96 

50 

3.582 

-22.9 

HOSLAK 

2738 

4 7 

11 

7.535 

240 

39 

43.542 

338.2 

SHEMAL 

2739 

52 

42 

55.761 

174 

6 

40.301 

43.4 

BfcLTSV 

2742 

39 

l 

39.845 

283 

10 

28.358 

8.0 

STNVIL 

2745 

33 

25 

31 .969 

2 69 

5 

10.453 

9.3 

CARGIL 

2809 

—46 

24 

43 .756 

168 

18 

13.551 

-0.3 

PARI BO 

2815 

5 

26 

53.104 

304 

47 

42.609 

-6.8 

MESHED 

2617 

36 

14 

26.218 

59 

37 

44. 326 

96 7.7 

FRTLMY 

2822 

12 

7 

53 .901 

15 

2 

6.953 

312.0 

NATLOP 

2837 

-5 

54 

57.951 

324 

49 

56.214 

30.0 

APL TWO 

2911 

39 

9 

48 .290 

283 

6 

14. 799 

112.2 

ETRPRE 

4050 

-25 

56 

37 .0 39 

28 

21 

29.238 

1573.8 

ETRMRT 

4082 

28 

25 

29.533 

279 

20 

6. 371 

-34.9 

NBER34 

4740 

32 

20 

53.441 

295 

20 

46.251 

-29.6 

NBER05 

4760 

32 

20 

53.662 

295 

20 

46.081 

-26.3 

NWAL 18 

4840 

37 

50 

29.394 

284 

30 

51.933 

-36 .4 

NWAL 13 

4860 

37 

51 

37.739 

284 

29 

24.704 

-37.6 

WC0R38 

494 6 

- 30 

49 

5 .327 

1 36 

50 

16.993 

68.3 

1UNDAK 

7034 

48 

l 

21 .267 

26? 

59 

20.298 

221.5 

1E0INB 

7036 

26 

22 

46.646 

261 

40 

7.997 

29.1 

ICOLBA 

7037 

38 

53 

36.121 

267 

47 

41.524 

236.0 

IBERMD 

7039 

32 

21 

49.581 

295 

20 

35.618 

-4.9 

I PUR 10 

7040 

18 

15 

28 .704 

294 

0 

24.056 

1 .9 

1GSECP 

7043 

39 

1 

15.530 

283 

10 

21.096 

13 .2 

IOENVR 

7045 

39 

38 

47.994 

255 

23 

39. 299 

1766.7 

GOULAS 

7050 

39 

l 

14.07? 

283 

10 

19.337 

12.9 

ROSLAS 

705 l 

35 

11 

47.130 

2 77 

7 

27.050 

848.3 
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TABLE 10. (continued) 


WALLAS 

7052 

w 

51 

35.886 

2 ttA 

29 

24. 034 

-28.2 

MGbL AS 

70S 3 


l 

15 .661 

283 

10 

19.604 

20.6 

CRML AS 

IQb** 

-24 

54 

15.095 

113 

A2 

58. 595 

20.8 

GM! SLS 

706 0 

1 3 

19 

33.753 

1 A A 

AA 

14. 244 

134.3 

STALAS 

706 3 

W 

l 

13.3A8 

28 3 

10 

20. 1 35 

14 .6 

ML030? 

706 5 

3 ) 

l 

1A .A('2 

283 

10 

19. 338 

14.1 

BD1LA1 

loti 

12 

21 

13.767 

295 

20 

38. 264 

-27.2 

GKKLAS 

iota 

21 

27 

37.756 

28P 

52 

5. 346 

-22.7 

1 JUM24 

7071 

2 7 

1 

1 3 .666 

279 

5 3 

13.276 

-18.5 

1 JUM40 

70 7? 

27 

1 

1A.262 

279 

53 

13.294 

-21.9 

l JtJPCI 

707 3 

?7 

1 

IA .508 

279 

53 

14. 199 

-30.1 

1 JU8C4 

70 7 A 

27 

l 

1 A . 7 1 A 

279 

53 

14.047 

-16 .8 

1SUDBR 

7075 

46 

27 

21.259 

2 79 

3 

11.097 

2 43.2 

1 JAMAC 

7076 

IB 

A 

3A.57A 

28 3 

11 

27.548 

427.9 

1GSFCN 

7077 

18 

59 

57.123 

28 3 

9 

38. 352 

12.5 

WALNOT 

101 a 

17 

si 

A7 .30A 

28A 

29 

28. 725 

-36.0 

1CARVN 

7070 

-24 

**4 

23 .16A 

113 

A3 

16. 916 

-1.3 

HAUL AS 

7800 

4 3 


56 . A 3 6 

5 

42 

45.037 

701 .2 

DAKLAS 

7820 

14 

46 

2.757 

3A? 

35 

28. 141 

42 .6 

GRASSE 

7842 

43 

45 

7.893 

6 

54 

10.914 

1306.3 

urglas 

700 1 

1? 

25 

2A .APR 

253 

26 

47.296 

1618.8 

OLILAS 

790? 

-2 r > 

57 

35 .8A1 

28 

14 

52. 921 

1566.2 

ARELAS 

790 7 

-16 

27 

56.733 

288 

30 

24.993 

2486 .0 

HfjPL AS 

7021 

11 

A 1 

3.187 

2 A9 

7 

19. 193 

234A.5 

NAT LAS 

79 2 9 

-8 

55 

A0.1A6 

32A 

50 

7.610 

32 .7 

grelas 

79 30 

3P 

A 

A2.17? 

2 3 

55 

57. 910 

506 .2 

OELFTH 

800 9 

s? 

0 

5.139 

A 

22 

16. 311 

28.6 

l IMWLD 

8010 

46 

52 

36.A22 

7 

27 

53.846 

931 .6 

MALVRN 

801 1 

52 

3 

35 .792 

358 

1 

54.174 

196.0 

HAUTtP 

8015 

43 

55 

57.017 

5 

42 

45. 175 

697 .4 

NICEFR 

8019 

4 3 

A3 

32 .A 32 

7 

17 

59.461 

420.0 

MUOONI 

80 30 

48 

A8 

21 .325 

2 

13 

46.859 

1 77 .8 

1 ORGAN 

900 l 

32 

25 

24.322 

253 

26 

49.217 

1623 .7 

1ULFAN 

900 2 

-25 

57 

36.021 

28 

1A 

52. 788 

1566 .9 

WGOMER 

900 3 

-31 

6 

2 .082 

i 36 

A 7 

3.407 

161 .2 

1SPAIN 

900 A 

36 

27 

A6.550 

353 

47 

37. 389 

67.6 

iTOKvn 

90 J 5 

15 

90 

22 .781 

1 39 

32 

16.897 

91.3 

1UATAL 

9006 

2s> 

21 

3A.5A1 

79 

27 

27.713 

1877.4 

1 UUI PA 

900 7 

-16 

27 

56.766 

288 

30 

24.889 

2489 .6 

1SMRAZ 

9008 

29 

38 

13.776 

52 

31 

11.325 

1587.8 

1 CUR AC 

9009 

12 

5 

25 .0 02 

291 

9 

44.866 

-17.1 

1 JUPTR 

9010 

27 

1 

13 .991 

279 

53 

13.788 

-16.3 

1 VI LUO 

901 1 

-31 

56 

3A . 7A0 

29A 

53 

36. 808 

632 .5 

1MAUI0 

9012 

20 

A? 

25.9AA 

203 

AA 

34. 363 

3048 .0 
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TABLE 10. (continuad) 


HOPMN 

9021 

31 

41 

2.961 

249 

7 

19. 151 

2347.8 

AUSPAK 

902 3 

-31 

23 

25.802 

1 36 

52 

43.906 

142.1 

DOOAIR 

9025 

36 

0 

19.853 

139 

11 

31.582 

896 .t 

DEZE1T 

902R 

8 

44 

50 .995 

38 

K 7 

33.831 

1912 .6 

NATALB 

9029 

-5 

55 

40 .443 

324 

50 

7.608 

39.1 

COMRIV 

90 n 

-45 

53 

12.521 

292 

23 

9. 807 

197.2 

JUPGEU 

9049 

27 

l 

14.039 

2 79 

5 3 

13.662 

-1 .0 

AGASSI 

9050 

42 

30 

21.398 

288 

26 

30.638 

147.5 

ATHENG 

905 L 

37 

58 

36.289 

23 

4 6 

40.689 

225.2 

MALVRN 

9000 

52 

8 

36 .017 

358 

1 

54.137 

166.1 

GREECE 

909 1 

38 

4 

44.203 

23 

55 

59.482 

503.9 

COLULK 

9424 

54 

44 

34.007 

249 

57 

23.449 

674.7 

EOWAF0 

9425 

34 

57 

50.636 

242 

5 

8.208 

754.7 

OSLOMR 

9426 

60 

12 

39.069 

10 

45 

3.968 

617 .2 

JflHNST 

9427 

16 

44 

39.021 

190 

29 

9.635 

26.6 

OLISAO 

9902 

-25 

51 

35 .894 

28 

14 

53.011 

1565.5 

ARESAG 

9907 

-16 

21 

56.753 

288 

30 

24.923 

2491 .6 

HOPSAO 

992 1 

31 

41 

3.189 

249 

7 

19.156 

2346.5 

MATS AO 

992 ) 

-5 

55 

40 .154 

324 

50 

7.562 

35.0 

GRESAO 

9940 

38 


42.412 

23 

55 

58.004 

503 .5 
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TABLE 11. A COMPARISON OF THE CALIBRATION AREA LASER STATIONS DETERMINED BY 
MARSH ET AL, 1977 WITH THOSE OF GEM 10 
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TABLE 12. THE MEAN EQUATORIAL RADIUS OF THE EARTH la e > 
DETERMINED FROM TRACKING STATION COORDINATES 



NUMBER 

ESTIMATED 

TRACKING 

OF 

a e 

SYSTEM 

STATIONS 

(meters) 

MOTS/SPEOPTS 

CAMERAS 

31 

6378140.8 

BAKER— NUNN 
CAMERAS 

27 

6378142.8 

LASERS 

16 

6378139.0 

DOPPLER 

23 

6378139.4 

ALL SYSTEMS 

114 

§378139.9 


U2 



where 


h - height of station above reference ellipsoid 
L - number of stations 

MSLH - mean sea level height of station from survey 
N - geoid height of station. 

A detailed gravimetric geoid mode! based upon GEM 10 was used for 
estimating the geoid heights. The short wavelength features of this geoid 
provided predominantly positive anomalies of 3 to 10 meters at a number of the 
stations, particularly those situated on Islands. The geoid was obtained privately 
from Marsh, and the method (which employed Stokes' function) is described in 
Marsh et al» 1976. 


3.5.2 a e Inferred from GEOS-3 Intensive Mode Altimetry 


GEOS-3 intensive mode altimetry was utilized for estimating the mean 
equatorial radius of the earth. The altimetry data set was selected for a 5° 
gridded distribution. These data were reduced in five day orbital arcs in which 
both laser and altimetry contributed to the determination of the orbit. GEM 10 
was used for the orbit and geoid computation. A single altimeter range bias was 
estimated from the altimeter residuals for each of these arcs. This altimeter 
bias contains all altimeter system biases along with the average error in the 
mean equatorial radius of the ellipsoid being used to compute the altimeter 
residuals. Martin (1977) has calibrated the intensive mode altimeter and finds it 
to measure short by 5.3 meters with a small uncertainty of 20 cm. Using this 
value for the system bias in the altimetry, the a e implied from ten five day arcs 
of altimetry is 6378141.0 m. These results are summarized in Table 13. A 
second important result to be noted in Table 13 is the exceptionally good fit to 
the altimeter data obtained using the GEM 10 geoid. GEM i0 did not use 
altimetry in its solution so this result can be viewed as a calibration of GEM 10. 
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TABLE 13. ESTIMATION OF THE MEAN EQ 
FROM GEOS— 3 INTENSIVE 


ARC EPOCH 

ARC LENGTH 
(DAYS) 

NO. OF 
ALT. C 

750616 

5 

2305 

750621 

5 

5454 

750527 

5 

2139 

750601 

5 

549 

750701 

5 

328 

750716 

5 

2251 

750730 

5 

3465 

750803 

5 

3555 

750815 

5 

2154 

750625 

5 

2229 


TOTAL/ 

AVERAGE 


24429 



a g IMPLIED 


ALTIMETER 
RMS OF FIT 
(METERS) 

BY RECOVERED 
ALT. BIAS 

6378000 + m 

2.27 

140.5 

2.57 

140.8 

2.22 

140.5 

2.67 

142.1 

1.64 

141.8 

2.78 

140.7 

2.04 

140.9 

2.50 

141.8 

2.55 

140.4 

2.63 

140.5 

2.46 

141.0 
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3.5.3 a Inferred From Mean Equatorial Gravity (g ) and GM 

G G 

A third method for estimating is based upon a new value of equatorial 
gravity g^, derived in GEM 10 from surface gravity data, and new value of GM. 
From the simple relation 


GM 



the variational relationship becomes 


a 

e 


or 


g " GM 




By using the old and new values of the reference parameters and bv 
removing the atmospheric mass ( = -87 x I0"S from the new satellite 

derived value of GM (398600.64 km /sec^), Ihe adjustment for a e is derived from 


the above equation. The result for a g and associated reference parameters^are 
given below in Table 14. The parameters refer to the old speed of light (C). 
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TABLE 14. GEODETIC REFERENCE PARAMETERS 


PARAMETER 

OLD 

NFW 

ACCURACY 

GM* 

398600.8 

398600.29 

+ .02 km^/sec' 

s e 

378031.0 

978031.5 

_r .5 mgal 

a e 

6378145 

6378139.3 

+ 1.5 m 


"Excludes the atmospheric mass and refers to c = 299792500 m/s. 


46 



3.5.4 Comparison of Fundamental GEM 10 Reference Parameters With 

Those Adopted by the IAG (1975) 

The GEM 10 reference parameters described in Section 3.5.3 are in 
remarkable agreement with the set established by a special study group of the 
IAG (Moritz, 1975). The GEM 10 values are adjusted to the IAG system by 
including in GM the atmospheric mass and the new speed of light (c = 299792458 
m/sec). The GEM 10 adopted a e is the composite value obtained from Sections 
3.5.1, 3.5.2, and 3.5.3 (a g - 6378140.). The GEM 10 values in the IAG system are 
compared to those adopted by the IAG in Table 15. The differences shown in 
Table 15 are very consistent with the uncertainties which have been stated for 
the parameters derived in the GEM 10 solution. 
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TABLE 15. COMPARISON OF THE IA6 1975 ANO GEM 10 
GEODETIC PARAMETERS 


PARAMETER 

IAG75 

GEM 10 

UNITS 

GM* 

398600.5 

398600.47 

km^/sec^ 

9 e 

978031.8 

978031.8 

mgal 

1/f 

298.257 

298.255 

— 

a e 

6378140 

6378140 

m 

c 

299792458 

299792458 

m/sec 


’Includes atmospheric mass (AGM = .35) and new speed of light (AGM = —.17). 
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SECTION 4, 


EVALUATION OF THE GRAVITY FIELD 




4. EVALUATION OF THE GRAVITY FIELD 


4.1 ERROR ESTIMATES OF THE POTENTIAL COEFFICIENTS 

It is alwa>s of interest to know the accuracy of a computed physical 
quantity as distinct from its formal precision measured by an experiment. In the 
case of comprehensive gravity model solutions, it has been shown by Lerch et al 
(1974) and Wagner (1976) that the formal uncertainties obtained from the 
solutions can be scaled to obtain reasonable estimates of the true uncertainties 
for the individual coefficients themselves. As indicated in equation 6 of Section 

3.2 a scale factor, f = \^T5^ was applied to the system of normal equations of 
GEM 9 and 10 in order to provide for realistic standard errors. We wish to test 
these error estimates here. Table 16 presents the coefficient errors for the GEM 
9 solution; the values in Table 16 represent the scaled error estimates 
(normalized) Tor the GEM 9 harmonics. Table 17 presents the estimated errors 
for the harmonics in the GEM 10 model. These error estimates were tested in 
three separate studies. We especially wished to confirm that a truly significant 
improvement has been obtained over previous GEM solutions for terms above 
degree 12 (which is indicated by significantly smaller uncertainties in GEM 9 and 
10 ). 


Rapp (1977) has estimated the terrestrial potential solely from surface 
data. Therefore, his model is completely independent of GEM 9 which was 
derived exclusively from satellite tracking data. Rapp f s model was used to 
calibrate the formal errors ascribed to the GEM 9 and 10 solutions. 

Figure 5 presents the estimated uncertainties from Rapp and GEM 9 
compared to the size of the coefficients from "Kaula's rule,” and those computed 
from these two solutions themselves. 
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TABLE 16. ESTIMATED ERRORS FOR OEM 9 (UNITS* 10~ V FOR FULLY NORMALIZED HARMONICS) 
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TABLE 17. EST MATED ERRORS FOR OEM 10 (UNITS - 10'* FOR FULLY NORMALIZED HARMONICS! 
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FIGURE 5. 

ROOT MEAN SQUARE POTENTIAL 



(GEM 9) 






Figure 6 is o calibration of the actual coefficient differences between 
the Rapp and GEM models compared to their estimated uncertainties. The low 
degree and order terms (/ < 12) in the Rapp model are not well determined from 
the gravimetry. The level of agreement between the uncertainties of the 
coefficients and the actual coefficient differences for the Rapp and GEM 9 
solutions (as exhibited in Figure 6) is remarkably good. It shows that the error 
estimates for the high degree terms (Table 16) are realistic. 

Surface gravity data (the 5° mean anomalies employed in GEM 10) were 
used as a second method to test the standard deviations of GEM 9. Commission 
errors (n^) of gravity anomaly due to errors in the GEM 9 model were derived 
from the gravity data based upon Kaula's statistics, (Kaula, 1966a). A scale 
factor f was computed to calibrate the standard deviations in GEM 9 as follows: 


Pk/m 1 * " 2( W 


where i = 978000 mgal. Results are given in Table 18 which are consistent for 
the various subsets of data and they verify the GEM 9 standard errors within a 
20% tolerance. The commission error of gravity anomaly based upon the GEM 9 
standard errors are plotted in Figure 7 as a function of the harmonics complete 
through degree t . 


A third approach using laser residuals was employed for testing the 
standard errors for the coefficients. The OR bitql AN alysis Program (ORAN, 
Martin, 1970) was used to integrate these coefficient errors as a gravity error 
model. The total estimated gravity error was propagated into simulated Grand 
Turk laser observations contained within a five day orbital reduction. High 
correlation in the errors of the zonal and resonance terms (m = 0, 13 and 14) 
required the elimination of their effects from the experiment. All other terms 
were included in the GEM 10 error model. 
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TABLE 18. CALIBRATION FACTOR (f) FOR GEM 9 STANDARD ERRORS BASED UPON 
COMMISSION ERRORS (o s > FROM 5° MEAN GRAVITY ANOMALIES 


NO. OF 5° 
MEAN 

ANOMALIES 

N ^ * 

a $ (mgal) 

f 

6 22 

25 

4.8 

1.1 

932 

20 

5.0 

1.1 

1109 

15 

5.2 

1.1 

1260 

10 

5.3 

1.2 

1404 

5 

5.6 

1.2 


*N is the number of 1° observed anomalies used in computing the 5° mean gravity anomaly. 
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FIGJJRE 7. COMMISSION ERROR OF GRAVITY ANOMALY 
BASED UPON GEM 9 STANDARD ERROR* 
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A five day orbit was computed using GEM 10 fitting the laser 
observations from August 4th to 9th, 1975. In this orbit, the laser data from 
Grand Turk (station No. 7068) was given zero weight and thereby did not 
contribute to the solution. The RMS (Root Mean Square of the residuals) fit to 
the Grand Turk observations, although unweighted in the solution, yielded an 
RMS of only 82 cm. In all, there were ten passes (2243 observations) of Grand 
Turk data. 

The estimated RMS predicted by ORAN for the GEM 10 gravity error 
contribution for all the Grand Turk measurement residuals was 78 cm. Gravity 
model error is the dominant error source in this test. The agreement between 
the ORAN simulation and the actual orbital fit to the Grand Turk data indicates 
that the standard errors for the coefficients are reasonable. An analysis similar 
to the above for estimating a gravity error model using ORAN is found in Martin 
and Roy, 1972. 


4.2 EVALUATION USING SURFACE GRAVIMETRY 

Surface gravity measurements are an important source of independent 
information for evaluating a global comprehensive gravity field. GEM 9 and 10 
hove been extensively "tudied using surface gravimetry. Figure 8 shows a 
c.mparison of recent GEM models with surface gravity. The GEM 9 field is in 
closer agreement with the independent surface gravity than any previous GEM 
satellite field. GEM 9 is in even closer agreement with this new surface 
'gravimetry than the GEM 6 (Lerch, et al, 1974) combination solution. GEM 10 
t also out-performs GEM 8. This is encouraging given the lower weight for the 
surface data (as discussed earlier) in GEM 10. Figure 9 compares recent surface 
gravity data sets with GEM 9. Quite clearly, the agreement between satellite 
and surface information is improving with time. 
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FIGURE 8. 


COMPARISON OF GEM MODELS WITH SURFACE GRAVITY 





Some recent combination solutions are: Goddard Space Flight Center - 
GEM 10, GEM 8, PGS 110 and GEM 8.1; the Smithsonian Astrophysical 
Observatory -SAO 4.3 (Gaposchkin, 1976); and GRIM 2 (Balmino, 1976). These are 
compared to the Rapp, 1977 surface gravimetry (Figure 10). GEM 8.1 is a repeat 
of the GEM 8 solution us'ng the new approach (described in Section 3.4) for 
combining the surface data but maintained the GEM d weight for the total data 
set. PGS NO is a repeat of the GEM 8 solution but complete to 30 x 30 instead of 
25 x 25. While GEM 10 performs very well, the relative weight of the surface 
data in GEM 8.1 was larger than GEM 10. Therefore, GEM 8.1 as would be 
expected agrees better with the surface data than does GEM 10. 

Table 19 shows the degree variances of the gravity anomalies from 
recent GEM solutions. The impact of the collocation (constraint) is noticeable in 
the loss of power in the high degree coefficients of GEM 9 and !0. The high 
degree coefficients in GEM 9 end 10 are somewhat smaller then their 
counterparts in recent GEM solutions. 


4.3 EVALUATION OF THE FREE AIR GRAVITY ANOMALIES DERIVED 

FROM GEM i AND 10 

A free air gravity anomaly map was computed from the complete GEM 9 
and 10 sets of coefficients. These maps are presented in Figures II end 12 
respectively. They are remarkably similar. Almost all gravity features are 
found in the same geographical location in these models, but hiere are occasional 
significant differences in the amplitudes for the indicated anomalies. Generally, 
when there is a significant difference in amplitude between the two fields, GEM 
10 shows anomalies with larger peak amplitudes. This is due to the surface 
gravity data providing greater definition of localized features. An example of 
this can be found over the Andes Mountains in South America. Both fields snow 
nearly identical placement for the anomah high in the Andes region, but in GEM 
10 the peuk is about 5 mgals larger. 
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FIGURE 10. 


COMPARISON OF COMBINATION MODELS WITH SUP FAC; GRAVITY DATA 


N 10 

1260 5° ANOMALIES (RAPP, 1976) 



GEOPOTENTIAL TRUNCATION DEGREE 
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TABLE 19. DEGREE VARIANCES OF GRAVITY ANOMALIES IN MGAL 2 


L 

,c /m + S 2 / m > 


DEGREE (/) 

GEM7 

GEM8 

GEM8.1 

GEM9 

GEM10 

3 

33.6 

33.7 

33.7 

33.5 

33.5 

4 

19.5 

19.6 

19.5 

19.5 

19.6 

5 

21.1 

21.1 

20.9 

20.7 

20.6 

6 

18.8 

19.1 

19.0 

18.9 

19.0 

7 

19.4 

18.3 

18.5 

19.3 

19.1 

8 

11.2 

10,2 

10.2 

11.6 

11.4 

9 

11.4 

11.2 

10.9 

11.4 

11.1 

10 

10.1 

9.8 

10.2 

10.0 

9.7 

11 

7.7 

7.3 

7.1 

6.7 

6.6 

12 

3.5 

3.2 

3.6 

3.6 

3.6 

13 

11.1 

6.5 

7.3 

6.5 

6.2 

14 

6.2 

3.3 

3.4 

4.0 

3.4 

15 

5.4 

4.5 

3.9 

3.2 

3.0 

16 

5.7 

3.5 

3.4 

2.3 

2.6 

17 

1.7 

6.9 

5.4 

2.0 

2.1 

18 

1.7 

5.0 

3.3 

3.3 

3.1 

19 

1.1 

9.4 

4.6 

2.9 

2.8 

20 

1.3 

8.4 

3.6 

2.2 

2.0 

21 

1.0 

5.9 

3.7 

1.1 

V8 

22 

0.5 

5.9 

3.4 

1.8 

1.7 

23 

0.2 

7.1 

3.1 

0.5 

0.5 

24 

0.2 

8.7 

3.8 

0.3 

0.4 

25 

0.6 

6.9 

3.3 

0.7 

0.7 

26 

0.7 

0.2 

0.4 

0.0 

0.0 

27 

1.0 

1.6 

0.9 

0.6 

0.6 

28 

^.6 

3.8 

3.5 

0.9 

0.9 

29 

2.1 

2.8 

2.0 

0.6 

0.6 

30 

0.0 

0.0 

0.9 

1.2 

1.2 
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There are basically two ways to independently assess the accuracy of a 
given gravity model. A direct comparison can be made between the satellite 
gravity models over areas where detailed surface gravimetry exists, (such as 
North America, Western Europe and Australia). A second comparison can be 
made between the geoid computed from the gravity model end the geoidal profile 
directly measured by sctellite altimeter experiments. This later approach is 
discussed in Section 4.4. The first approach is discussed below. 

Since the low degree portion of the gravity fields are generally 
recognized as being accurately determined from satellite observations, we have 
concentrated our comparisons on the higher order terms in the model. Figure 13 
presents the geographical distribution of the surface data in Rapp's (1977) 
potential model computed solely from surface data. Rapp's data set of 1° x 1° 
free air anomalies cover approximately 68% of the earth's surface. Almost two- 
thirds of the measurements are in the northern hemispere, however. 

Figure 14 presents a map of the Rapp free air anomalies computed for 
coefficients of degree 13 to 22 from his model. The contour interval is 4 mgals. 
The darker areas are those where the free air anomalies are less than -4 mgals. 
The lighter shaded regions are areas with small gravity signal at this wavelength 
being from -4 to +4 mgals. The white areas locate positive anomaly features 
being greater than 4 mgals. The half wavelength for this portion of the gravity 
field ranges from 1500 to 900 km. 

Figure 15 is a similar free air anomaly map from GEM 9 for coefficients 
of degree 13 to 22. The GEM 9 model is completely independent from the model 
computed by Rapp since it uses only satellite tracking data. Figure 16 overlays 
the boundaries of the Rapp inferred anomalies onto the GEM 9 anomaly map. 
The agreement in terms of the geographical location of the anomalies is striking. 
Those areas which have good gravimetry show excellent agreement between the 
Rapp model and GEM 9. In almost all instances, the discrepancies in this 
comparison occur in those regions where Rapp does not have data (e.g., the 
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FIGURE 15. FREE AIR GRAVITY ANOMALIES FOR THE GEM 9 
GRAVITY MODEL FOR COEFFICIENTS OF DEGREE 13 TO 22 
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southern oceans). This intercomparison demonstrates that the satellite derived 
gravity models are becoming increasingly more accurate in their ability to 
resolve relatively short wavelength gravitational features. The comparisons with 
altimetry presented irf Section 4 .4 confirm this conclusion. 

The relatively high degree portion ( l = 13 to 22) of the satellite derived 
gravity field is of geophysical interest. Therefore, we have prepared a map 
(Figure 17) of the estimated gravity anomalies of the upper mantle derived from 
GEM 10 with the crustal features removed (using the isostatic model of Khan 
(1973)). This map shows the estimated mantle gravity features of half 
wavelengths ranging from 1500 to 900 km. To facilitate on analysis of 
convective processes, we have indicated the tectonic plate boundries obtained 
from Chappie and Tullis, 1977. 


4.4 EVALUATION OF GEM 9 AND 10 USING ALTIMETER DATA 

4.4. 1 Evaluation of GEM-9 and 10 Using the "Round the World 11 Data Taker. 

from Skylab 

The SKYLAB- 1 93 radar altimeter was operated nearly continuously 
around the world on January 31, 1974. This direct measurement of the sea 
surface topography provided for the first time an independent basis for the 
evaluation of a global geoid computed from satellite derived gravity models. The 
models considered were the Goddard Space Flight Center GEM (1-10) models; the 
Smithsonian Astrophysical Observatory SAO 4.3 model, and GRIM 2. This data 
has previously been used by Marsh et. al., 1975 for gravity model evaluation. The 
results obtained in our analysis differ somewhat from those of Marsh. A time tag 
error was discovered in the application of the SKYLAB Airlock Module Time. 
This error bias has been corrected in our tests. The "round the world" data 
consisted of 396 six second smoothed altimeter ranges which encircled the world. 
The RMS of fit to this data is shown in Figure 18. The 3.16 and 3.01 meter 
residual RMS from GEM 9 and 10, respectively, is quite satisfying. Contained 
within these residuals are: 
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GEOlD MT —METERS 


FIGURE 18. SKYLAB "ROUND THE WORLD" ALTIMETER 
RESIDUAL RMS BEFORE CORRECTION FOR A TIME TAG ERROR 



SKYLAB "ROUND THE WORLD" DATA TAKE (RMS - 8m) 
31 JANUARY 1974 


ALTIMETER RESIDUALS AFTER TIME TAG CORRECTION 


MODEL 

SURFACE GRAVITY ONLY 

SA04.3 

GRIM2 

GEM 1 
GEM 2 
GEM 3 
GEM 4 
GEM 5 
GEM 6 
GEM 7 
GEM 8 
GEM 9 
GEM 10 


RMS OF 396 ALTIMETER 
OBS. IN METERS 

6.25 

6.21 

5.70 

3.74 

3.91 

4.08 

5.13 

3.89 

4.47 

3.28 

4.57 

3.16 

3.01 
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• commission error in the computntion of the sea .'face from 
the GEM 9 and 10 coefficients themselves, 

• omission error in the same computation from the models c ;e 
to their truncation, 

• altimeter noise which for SKYLAB was assessed to be I to 2 
meters, and 

• orbital error in the radial positioning of SKYLAB. 

The truncation error by itself is estimated to be nearly 2.5 meters. 


4.4.2 Evaluation of GEM 10 using GEOS-3 Intensive Mode Altimeter Ogta 

The GEOS-3 altimeter was c^erally operated over specific geographic 
areas dicing a specified period of a few weeks. These areas were varied over 
time so that a global data set could be compiled from the total complement 
of acquired GEOS-3 altimeter passes. This type of data accumulation does 
not lend itself to the global calibration of a gravity model. The time requiring 
precision orbit determination with all the data would be about a year. However, 
during February and March of 1976, the altimeter was operated in a continuous 
fashion and a reasonable, although not completely global, distribution of altimetry 
is available. 

A test to independently assess the accuracy of the geoid from some 
recent gravity models was designed. Two five-day orbital arcs were reduced 
during this concentrated tracking period. The first extended from February 
29 to March 4, 1976 while the second was from March 10 to 15, 1976. Each 
orbit determination made use of all the laser and intensive mode altimeter 
data available during these intervals. Figure 19 shows the distribution of the 
altimeter passes which were employed in the two solutions. 
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The five day orbits were determined using the complete GEM 10, SAG 

4.3 and GEM 8 gravity models. The reference geoid used in computing the 

altimetei measurement residuals came from these models respectively. Only 

intensive mode altimetry was used; a sirgle bias representing the mismodeling of 

a and instrument bias was solved for in each of the 5 day arcs, 
e ' 

Table 20 summarizes the RMS of fit to the 10750 altimeter observations 
contained within the two five day arcs. The GEM 10 results were excellent; they 
showed an even closer agreement between the GEM 10 geoid and the altimeter 
profiles than was seen in the SKYLAB comparisons (Section 4.4.1). The GEM 10 
results are completely consistent with those presented in Section 3.5.2 for 10 
five day arcs having less globally distributed data. The way this test was 
performed makes it difficult to attribute the poorer resuits obtained from SAO 
4.3 and GEM 8 to geoid error. Orbital error is also contained within the residuals 
from the respective models and probably contributes a sizable amount to the 
total residual RMS obtained. 

On the other hand, analysis of 42 short arcs (10-20 min./arc) data in 
which the orbital errors were removed empirically show significant improvement 
in the geoid from GEM 9 and 10 over GEM 7 and 8. Residuals of altimeter 
derived sea surface with geoid heights from these fields in these globally 
distributed arcs were 3.30 m for GEM 8, 2.85 m for GEM 7, 2.66 m for GEM 9 
and 2.52 m for GEM 10. The improvement of GEM 9 over GEM 8 is especially 
gratifying since GEM 9 is a smaller field without the benefit of surface data. 

In a second case, the GEM 10 field was truncated at twelfth degree and 
order for the computations of the geoid, while the orbit determined previously 
from the full GEM 10 field was retained. This variation of the test wos made to 
assess what degradation, if any, would result in eliminating the contribution of 
the high degree and order terms to the GEM 10 geoid. The altimeter residual 
RMS (Table 12) increased by almost 1.5 meters when these higher degree and 
order terms were eliminated. The degradation due to truncation of GEM 10 to 
12x12 can be estimated by: 
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TABLE 20. 

GEOS— 3 INTENSIVE MODE ALTIMETER 
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L 


1/2 



wltere 

L is the loss of accuracy due tc truncation 

Tj is the total RMS from the truncated solutions (combined), and 

Tp is the total RMS from f he solutions using the full model 
(combined). 

Combining the results from the two cases, the estimated loss of accuracy 
due to the truncation of GEM 10 to 12 x 12 is 3.22 m. This is compelling 
evidence that the higher degree and order coefficients in GEM 10 contribute 
accurate information to the computation of a global geoid. This strengtnens the 
conclusions made in Section 4.3. The worldwide contribution of the terms of 
degree 13 to 22 in GEM 10 is at least 2.5 meters (rms). This would further 
indicate that this portion of the GEM 10 model is highly accurate. 


4.5 EVALUATION OF GEM 9 AND 10 USING ATS-6/GEOS-3 DOPPLER 
EXCHANGE DATA 

GEOS-3 and ATS-6 performed a four-way doppler exchange experiment 
(Satellite to Satellite Experiment: SSE). The SSE data were not included in GEM 

9 and 10. A two revolution orbit of GEOS-3 (revolutions 245 ar.d 246) was 
reduced which had particularly strong laser ground tracking and two consecutive 
45 minute SSE tracks* Figure 20 summarizes these results. This test was 
designed to evaluate the high frequency portion of the geoootentia! model. GEM 

10 fit 4 he dote particularly well. The randomness (RND) of the SSE residuals is 
also listed. GEM 10 again was the superior solution though a small signal still 
remains in these residuals. 
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RESIDUALS (CM/SECI RESIDUALS (CM/SEC) 


FIGURE 20 
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RESIDUALS 

normal j 

MODEL 

CM SEC 

- DEVIATE ' 

GEM9 

096 

19 7 

GEM10 

089 

16 0 

GEM7 

131 

54 8 

GEMS 

108 

29 3 

SAO 4 3 

233 

221 8 

GRIM2 

197 

159 9 
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4.6 


EVALUATION OF 13th ORDER HARMONICS USING RESONANT 
SATELLITE ORBIT > 


Klosko and Wagner (1975) used over 130 constraint equations developed 
from the analysis of deep resonance orbital passages, new shallow resonance 
harmonic determinations, and the frequency decomposition of existing satellite 
geopotential models to obtain improved values for the 13th order tesserol 
harmonics. In all, thirteen satellite orbits having inclinations from 2 1‘° to 
retrograde were evaluated for this solution. The estimated harmonics were 
complete to the 32nd degree. 

The 13th order coefficients obtained from this resonance solution are 
compared with GEM 7 and GEM 9 (Figure 21 and Table 2IL GEM 9 is in closer 
agreement with the resonance information than was GEM 7. Term C28,I3 seems 
to show large variation* f rom solution to solution. When this term is removed 
from the comparison, the GEM 9 field has less than one half the RMS for 
coefficient differences when compared to the resonance solution than had GEM 
7. This result is all the more surprising because the shallow resonance 
information GEM 7 is a significant componen! of the Klosko and Wagner 

solution. 
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TABLE 21. 

COMPARISON OP 13TH ORDER COEFFICIENTS FROM GEM MODELS WITH 
THOSE DERIVED FROM RESONANCE ANALYSIS (NORMALIZED VALUE x 10 9 ) 


r 

t — 

m = 13 

(D 

KLOSKO, 
WAGNER 1975 
NORM. VALUE 

0 

(2) 

GEM7 

(3) 

GEM9 

1 - 2 

1 - 3 

C 13 

-62.6 

( .9) 

-60.9 

-60.1 

- 1.7 

- 2.5 

S 13 

68.6 

( .9) 

67.6 

69.8 

1.0 

- 1.2 

C 14 

29.5 

(2.1) 

256 

28.0 

39 

1 5 

S 14 

46.9 

(1.1) 

43.9 

422 

3.0 

4.7 

C 15 

—24.6 

(1.2) 

-25.3 

-22.8 

0.7 

- 1.8 

S 15 

- 4.9 

(1.0) 

- 3.1 

- 2.2 

- 1.8 

- 2.7 

C 16 

16.7 

(1.8) 

8.4 

12.2 

8.3 

4.5 

S 16 

1.7 

(1.3) 

1.6 

- 7.5 

0.1 

92 

C 17 

14.4 

(1.9) 

16.8 

14.7 

- 2.4 

- 0.3 

S 17 

18.7 

(1.8) 

24.9 

19.2 

- 6.2 

- 0.5 

C 18 

- 5.8 

(2.8) 

-19.2 

-12.0 

13.4 

62 

S 18 

-32.7 

(1.7) 

-35.8 

-37.4 

3.1 

4.7 

C 19 

-12.7 

(3.8) 

- 9.4 

-12.4 


- 0.3 

S 19 

-30.6 

(3.7) 

-20.2 

-30.7 

-10.4 

0.1 

C 20 

21.9 

(3.2) 

2.9 

23.2 

19.0 

- 1.3 

S 20 

5.0 

(2.3) 

8.2 

3.9 

- 3,2 

1.1 

C 21 

-15.0 

(3.5) 

-21.9 

-16.4 

6.9 

1.4 

S 21 

9.7 

(3.7) 

16.2 

13.9 

- 6.5 

- 4 2 

C 22 

-17.5 

(3.0) 

-24.3 

-30.0 

6.8 

12.5 

S 22 

12.0 

(2.5) 

19.3 

7 7 

- 7 3 

4.3 

C 23 

3.8 

(4.0) 

- 7.0 

- 2 2 

10.8 

6.0 

S 23 

- 4.3 

(4.1) 

4.9 

- 1.6 

9 2 

- 2.7 

C 24 

1.4 

(4.9) 

-14.9 

6.1 

163 

- 4.7 

S 24 

- 2.1 

(3.7) 

3.8 

- 7 3 

- 5.9 

5.2 

C 25 

10.1 

( 7.0) 

25.1 

152 

-15.0 

- 5 1 

S 25 

-11.6 

(6.3) 

11.1 

- 8 8 

-22 7 

- 2 8 
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TABLE 21. (continued) 

COMPARISON OF 13TH ORDER COEFFICIENTS FROM GEM MODELS WITH 
THOSE DERIVED FROM RESONANCE ANALYSIS (NORMALIZED VALUE x 10 9 ) (cont.) 

( 1 ) 


c - 

m = 13 

KLOSKO r 
WAGNER 1975 
NORM. VALUE 

0 

(2) 

GEM7 

(3) 

GEM9 

1 - 2 

1 - 3 

C 26 

1.3 

(7.5) 

5.1 

- 3.2 

- 3.8 

4.5 

S 26 

0.7 

(5.5) 

4.0 

- 8.2 

- 3.3 

8.9 

C 27 

- 7.0 

(7.7) 

16.4 

- 7.7 

-23.4 

07 

S 27 

- 6.8 

(6.9) 

21.2 

-10.7 

-28.0 

3.9 

C 28 

-16.5 

(8.0) 

-55.9 

20.6 

39.4 

-37.1 

S 28 

0.4 

(5.2) 

9.2 

8.8 

- 8.8 

- 8.4 

C 29 

-17.4 

(5.6) 

- 6.5 

-10.8 

-10.9 

- 6.6 

S 29 

- 9.5 

(5.3) 

- 7.4 

- 9.0 

- 2.1 

- 0.5 

C 30 

9.9 

(8.1) 

— 

— 



S 30 

4.7 

(4.0) 

— 

— 





RMS OF RESIDUALS 

w/C 28 : 

12.56 

7.93 





w/out 

Co®: 

10.75 

4.80 
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SECTION 5. 


EVALUATION Or GEM 9 AND 10 FOR ORBIT 
DETERMINATION ACCURACY 




5. EVALUATION OF GEM 9 AND 10 FOR ORBIT 
DETERMINATION ACCURACY 


The ability to model accurately the gravitational forces on near earth 
satellites is one of the most important applications for improved geopotential 
models. GEM 9 and 10 have undergone extensive testing in this regard. 

The GEOS-3 orbital accuracies were of paramount concern given the 
demands of altimeter support. But also of concern was the quality of the 
computed orbits for Beacon Explorer-C (BE-C) and LANDSAT. LANDSAT data 
was used for the first time in GEM 9 and 10. 

BE-C is used extensively in the Laser Polar Motion and San Andreas 
Fault Experiments. As such, it has been extensively tracked by various laser 
systems. Table 22 presents the results obtained using laser data from BE-C. 
This laser data is not in the GEM 9 or 10 solutions, though other range data are 
used from the same stations to BE-C. Two station configurations -those on the 
East coast of the United States and those on the West coast - were tested. GEM 
9 and 10 show considerable reduction in the overall fit to this laser data when 
compared to other available models. In the case of Storiette, Marsh and 
Williamson (1976) have extensively analyzed the orbital accuracies obtained from 
some preliminary GEM models. 

The GE'OS-3 spacecraft is in a nearly circular orbit at an altitude of 

approximately 840 km. The spacecraft is not extremely dense and has an 

area/mass ratio of 1 .4365m^/345.909 kg (.004). At this altitude, the estimates of 

2 

the atmospheric drag perturbations on GEOS-3 range from I2m/day to 

2 

20m/day (along track) when using the Jacchia (1971) Density Model. A drag 
perturbation of this magnitude requires extremely refined modeling to avoid 
prohibitively large orbital positioning errors. We account for the drag on GEOS-3 
in a variety of ways depending on the length of the orbit to be determined; 
briefly: 


PRECEDING PAGE BLANK NOT FILMS) 
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TABLE 22. 

BEACON EXPLORER-C (BE-C) 

LASER RESIOUAL RMS FROM THE SAFE EXPERIMENT 


MODEL 

EAST COAST 
OBSERVATIONS 
(3 SIX-HOUR ARCS) 
LASER RMS IN CM 

WEST COAST 
OBSERVATIONS 
(2 FOUR-HOUR ARCS) 
LASER RMS IN CM 

GEM 7 

54 

50 

GEM 8 

126 

39 

GEM 9 

18 

23 

GEM 10 

29 

18 

SA0 4.3 

280 

154 

GRIM2 

756 

269 
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• When the orbit is less than 12 hours in length, a ballistic 
coefficient (Cq) is modeled at a fixed value of 2.5. Any small 
residual error is easily absorbed .n the epoch parameters of the 
orbit. 

• When the orbit is longer than 12 hours but shorter than 36 hours, 
Cq is allowed to adjust. 

• Lastly, when the orbit is longer than 36 hours, a time varying as 
well as a constant Cq are adjusted to the data. 

The GEODYN Program (T. Martin, 1972) is used for the orbital 
reductions. GEODYN uses Cowell type numerical integration techniques. For 
GEOS-3 orbital reductions, luni-solar gravitational perturbations, solar radiation 
pressure, BIH polar motion and UT I data and atmospheric drag using the Jacchia 
1971 Density Model are modeled. We also model solid earth tides ((<£ = .29) and 
the ocean tides using the diurnal lunar model of Hendershott (1970). 

The orbits calculated for GEOS-3 were thoroughly tested. The radial 
accuracy of GEM 9 and 10 has been evaluated using intersecting GEOS-3 
altimetry passes from independent and widely separated orbits (in time). The 
crossover points were differenced to estimate the radial error in the GECS-3 
orbits computed from GEM 9 and 10. 

The altimetry residual for the Kth revolution is given by 

res K = a K ' (r K " 9 K " V 

where 

a = altimeter range 

r = satellite height above the reference ellipsoid 

g = geoid height 

t = tide height 
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V 


Since the residuals are differenced over the same location, the geoid 
height cancels. Tides were modeled using Hendershott (1970) and small errors are 
present. Ignoring the tides (t), 'he difference of the satellite altimetry residuals 
for the Kth and Jth po ' at intersection K, J is 

Ares K,J H (a K " a J } - (r K * r J ) * 

If the altimetry is assumed to be noiseless, and having a constant bias, then the 
altimeter crossover residual difference is a measure of radial orbital error. 

Four one-day arcs spanning the altimeter measurements were computed 
from laser range data; the altimetry was not used in the determination of the 
orbit. There were II intersections in t e "Atlantic Ocean region and 28 
intersections southwest of Australia as illustrated in Figure 22. T he altimetry 
intersections in the Atlantic region involved at least one pass of altimetry in the 
global mode which Martin (1977) has shown to have varying off-nodir biases of 
from Mo 3 meters. The Australia intersections were computed from altimetry 
which was all in the intensive mode. The intensive mode data has a known bias 
of a constant (-5.30 + .21 meters) but is not noticeably affected by pointing 
errors. These data therefore should yield superior crossover results. 

Table 23 lists the crossover results obtained from GEM 9 and 10 a : ong 
with other representative fields. The RMS is given separately for intensive 
mode, global mode and the total set of 39 intersections. It is readily seen that 
the intensive mode is much more arcurcfe than the global mode even though the 
Australia intersections are further from tracking stations. The radial error for 
GEM 9 and 10 in these tests appears to be less than I meter. 
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TABLE 23. 

RMS OF ALTIMETER CROSS OVER RESIDUALS (METERS) 
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FIGURE 23. 

DATA DISTRIBUTION OF FOUR ARCS USED IN GEM10 LONG ARC ANALYSIS 
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A second altimeter cross over test was comoleted using the GEM 10 
gravity model. In this study four long orbits (three 5 day and one 4 day arc 
lengths) were determined from available laser data. Figure 23 gives the details 
of the laser tracking data used in these trajectories. Arcs one end tv*o used 
tracking almost exclusively from the NASA lasers in the GEOS-3 Calibration 
Area. Arc two was especially weak since the laser data at Patrick AFB 
(RAMLAS) had timing problems; RAMLAS timing biases had to be estimated 
from this data simultcneously with the orbit thereby further reducing the 
strength of the solution. Arcs three and four had a good distribution of NASA 
and SAO laser data. The altimeter data (intensive mode only) was not used in the 
orbit determinations. 

Figure 24 shows the location of the 127 altimeter cross over points 
obtained by intercomparing all four of these arcs. The cross over distribution is 
still unbalanced, but it is nearly global for sampling different parts of the orbit 
especially those parts away from the tracking stations. 

Figure 25 presents a histogram of the GEM 10 altimeter cross aver 
residuals. Arc two has been segregated by itself and the results from this a*"c do 
show an anticipated degradation in radial orbital accuracy. The 80 intersections 
which do not involve data from arc two have a residual RMS of 1.31 meters. The 
total RMS (including arc two) for 127 cross overs is 1.60 meters. These result;, 
reflect numerous errors besides radial orbital errors. 

An error budget for the crossover results is estimated by: 


R 2 = 2(E 2 ) ♦ G 2 ♦ 2(T 2 ) ♦ 2(^ 2 ) 
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FIGURE 24. LOCATION OF 127 ALTIMETER 
INTERSECTIONS FROM THE GEM 10 LONG 
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FIGURE 25. HISTOGRAM OF GEM 10 LONG ARC ALTIMETER CROSS OVER TEST 



CONTAIN FEB. 1976 ARC 
IN CROSS OVER (WITH 
RAMLAS TRACKING) 



□ DOES NOT CONTAIN 
FEB. 1976 ARC IN 
CROSSOVER 


TOTAL RMS FOR 127 INTERSECTIONS 
1.60 m 

RMS FOR 80 INTERSECTIONS WHICH 
DO NOT CONTAIN FEB 1976 
ARC = 1.31 m 


iiimiiiiHU'im X 


Oto 1 


1to2 2 to 3 3 to 4 4 to 5 5 to 10 10 and over 


ALTIMETER CROSS OVER RESIDUAL IN METERS 


96 



where 


R is the total residual altimeter cross over RMS; (l.3lm) 

G is the geoid height mismatch. The crossover data were compiled 

by hand and the altimetry was not interpolated to obtain a value 
at the precise intersection point. Rather, the closest points in the 
respective passes were used and these can be specially separated 
by as much as 20 km. An estimate of this mismatch is .3 m on 
average. 

T is the ocean tidal error. This has been estimated to be .3 m in 
each arc, on average from the Hendershott Model. 

£ is the altimeter noise. In these tests, we used the uneditted major 
frame averages made available from Wallops Space Flight Center. 
Our noise estimate is .3 meters in each arc, and 

E is the orbital error in each arc. 

When this equation is solved, the estimated orbital error is .80 m from 
the three long arcs. With all tour arcs, the estimated radial orbit error is 1.03 m. 

The GEOS- 3 orbital accuracies from GEM 9 and 10 have also been 
extensively tested on 3 revolution, I day and 5 day arcs estimated from laser 
range data. Appendix I presents these results. 

Tests I through 9 show various methods employed in determining the 
accuracy of the fields. The basic approach has been to intercompare two 
different orbit trajectories. For example, an orbit is determined over a period of 
time. A shorter arc length within the first is selected and its trajectory is 
determined. We then compare the two solutions in their radial, crosstrack and 
along track component differences over their common interval. In this way we 
can evaluate the accuracy of the field in all three components. 
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In tests I through 5, five day orbital arcs are compared with one and two 
day arcs which all contain subsets of the same data. Test 6 navigates a station 
height by using separate passes of data not in the orbital solution. Since we 
allow the station to adjust only in height, it is a good evaluation of the radial 
accuracy of the field for high elevation passes. We chose one northern and two 
southern hemisphere stations to insure that we had a good global sample. 

In tests 7 and 8 we did not include the data from the short arc in the 
longer one. This ieft a gap of from one day to 32 hours in the longer arc. Test 9 
evaluates the RMS of fit for a 5 day orbital arc determined from laser data. All 
tests are described in detail on their individual summaries. 
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SECTION 6 . 


SUMMARY 




PRECEDING PAGE BLANK NOT FIT. MED 

b. SUMMARY 


The major objectives of GEM 9 and 10 were achieved. GEOS-3 orbital 
accuracies from these models are about I m in their radial components for 5 day 
arc lengths. The new GEM 9 and 10 models yield significantly improved results 
when compared to the surface gravimetry, SKYLAR and GEOS-3 altimetry and 
highly accurate BE-C laser ranges than do previous GEM solutions. We believe 
that a genuine improvement has been achieved for the global representation of 
the terrestrial potential. 

Additionally, a new value of GM has been determined dynamically from 
laser tracking. A consistent value of the mean equatorial radius of the earth was 
obtained from the estimated tracking station coordinates, the GEOS-3 altimeter 
data and the • nplied value of g e . The average value of a e was found to be nearly 
constant among the different techniques used to estimate this parameter. The 
set of recommended or adopted physical constants from this work are: 

• GM = 398600,64 km^/sec 2 

• a = 6 n7 8l40m 

e 

• c r 299792.5 km/sec 

• f = 1/298.255 

The accuracies of the geopotential coefficients have been estimated end 
imply commission errors in geoid height of 1.9 m and 1.5 m (global RMS values) 
respectively for GEM 9 and 10. This error estimate was obtained from 
independent calibrations with the surface gravimetry, the GEOS-3 altimetry and 
an error propagation using a gravity model error model derived from these 
estimates. 
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APPENDIX I 


GEOS- 3 ORBITAL TESTS 




APPENDIX 1. 

GEOS 3 OR31TAL TESTS 
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LASER DATA SET FOR GEOS-III 
(APRIL THROUGH AUGUST 1975 ) 
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TEST * DESCRIPTION: A five day arc is determined from laser data. Within the same data span, five 
one day arcs are also determined from the laser data. These orbital trajectories are then differenced 
every minute over their common time interval, ^he radial (R), cross track (C) and along track (A) 
position differences are statistically evaluated as an RMS difference in each of these ballistic components. 
These RMS differences for various test gravity models for each comparison (i.e.. the 5 day arc versus each 
one day arc) are presented. The 5 day arc selected (May 18 to 23, 1975) was viewed by us as being very 
weak given the limited amount of tracking data available. Some of the one day arcs are also ve*y weak. 


ORBITAL COMPARISON: RV 
5/18 to 23/75 

5° ARC VERSUS EACH ONE DAY ARC 


GRAV. 

MODEL 


5 D v.1? 

5°v. 1° 

5° v. 

i° 

'3 

5° v. 1? 

5°v. 


R 

4.23 

6.83 

5.23 


„.40 

11.32 

GEM7 

C 

15.91 

14.04 

3.98 


7 29 

11.71 


A 

13.44 

16.18 

42.84 


23.45 

29.32 


R 

2.65 

1.24 

0.88 


0.21 

1 84 

GEM9 

C 

9.04 

0.42 

1.52 


1.43 

1.83 


A 

6.44 

3.92 

8 79 


0.84 

5.37 


R 

2.64 

0.69 

0 79 


0.11 

1.74 

GEM10 

C 

8.92 

1.03 

148 


1.32 

1 77 


A 

644 

2.15 

8.35 


0.72 

4.26 

NO. OF 
PASSES 


2 

3 

4 


5 

3 


IN EACH 
1 D ARC 


D ^ DAYS 
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TEST 2 DESCRIPTION: A five day arc is determined from laser data. Within the same data span, 
a two day arc. and two one day arcs are also determined from the laser data. These orbital trajec- 
tories are then differenced every minute over their common time interval . The radial (R), cross 
track (C) and along track (A) position differences are statistically evaluated as an RMS difference m 
each of these ballistic components. These RMS differences for various test gravity models for each 
comparison (i.e., 5 day versus 2 day arcs. 5 day versus 1 day arc No. 1 , S day versus 1 dsy arc No. 2, 
etc), are presented. The 5 day arc selected (May 18 to 23, 1975) was viewed by us as being very 
weak given the limited amount of tracking data available. The second one day arc <5/22 to 23 > had 
only 3 passes of data, none of whicn was past the 11th hour on this day while this orbit was 
differenced for a full 24 hour interval. 


ORBITAL COMPARISON RV 
5/18 TO 23 75 

5° ARC; 2° ARC *5 21 - 23). *5 21 i; 5 22 • 


GRAV. 

MODEL 


5 D v.2 D 

5°v.1? 

5°v.1° 

2°v.1? 

2°v 1 


R 

4.70 

840 

11.32 

3.93 

11.71 

GEM7 

C 

6 63 

7.29 

11.71 

1.43 

15 05 


A 

23 58 

23.44 

29 32 

10.15 

39.94 


R 

0.79 

2.24 

4.96 

1.62 

5.15 

GEM8 

C 

4.68 

2.26 

8.19 

3.47 

3,79 


A 

22.64 

7,92 

61.51 

4.16 

31.32 


R 

4.79 

5.36 

28.97 

0.59 

32.18 

GRIM2 

C 

8 29 

15.41 

40.35 

23.10 

45.57 


A 

10.12 

12.71 

>50 

5.26 

>50 


R 

5.08 

7.67 

44.04 

4.89 

49.02 

SA04.3 

C 

8.19 

14.47 

63.69 

19 53 

71.22 


A 

28 74 

19.35 

>50 

12.95 

>50 


R 

0.34 

0.21 

1 84 

0.54 

1.63 

GEM9 

C 

0.24 

1.43 

1.83 

1.19 

1.84 


A 

2.26 

0.84 

536 

1.39 

7.55 


R 

0.42 

0.11 

1.74 

0.53 

1.41 

GEM10 

C 

0 29 

1.32 

1.77 

1.03 

1.64 


A 

2.36 

0.74 

4 25 

1.29 

6.04 


D 

2 


COMMENT 


*12 has only 3 passes of data to 1 1 h and 13 h predict. 
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TEST 3 DESCRIPTION: A five day arc (August 2-7. 1975) was determined from laser data. 
A one day arc (August 3—4. 1975) was determined from the same laser data. These orbital 
trajectories are then differenced for every minute over their common one day period and an 
RMS is computed for the difference in each ballistic component — (R) radial. <C) cross track, 
(A) atom? track. 1 his arc was viewed by us as being exceptionally well tracked from the laser 
system and should yield strong orbit determination possibilities. 


ORBITAL COMPARISON: RV 
8/2 TO 7/75 
5° ARC; 1° 8/3 - 4 


GRAV. 

MODEL 


5 D v.1° 


GEM7 


R 6.41 

C 305 

A 31.67 


GRIM2 


R 10.15 

C 2.26 

A 38 23 


SA04.3 


R 7.56 

C 1686 

A 36.34 


GEM9 


R 053 

C 0.43 

A 2 23 


GFM10 


R 0.47 

C 044 

A 2.03 


•1 


D 


has 11 passes of User data. 


10 ° 



TEST 4 DESCRIPTION: A five day arc (October 27 to November 1, 1975) was determined 
from laser data. A two day arc (October 29 to 31) and a one day arc (October 29) was aho 
determined from the same laser data. These orbital trajectories are then differenced for every 
minute of their common periods and an RMS is computed for the difference in each ballistic 
component. 


ORBITAL COMPARISON: RV 
10/27 TO 31/75 
5°: 2° (29 - 311; 1° (29 - 30) 


GRAV. 

MODEL 


5^.2° 

5° v. 1 D 

2° v. 


K 

7.40 

7.31 


GR IM2 

C 

23.53 

20.92 



A 

44.32 

49.85 




R 

4.38 

8.45 

5.03 

3A04.3 

C 

40 20 

42.37 

5.15 


A 

72.12 

46^9 

31.81 



R 

1.99 

3.59 

3.47 

GEM7 

C 

5.64 

647 

2.43 


A 

$.88 

26.42 

20.72 



R 

7.09 

7.49 

GEMS 

C 

4.€4 

4.87 


A 

**0.17 

42.43 



B 

046 

0.69 

0.27 

GEM9 

C 

1.25 

1.47 

0.23 


A 

2.33 

2.33 

2.87 



R 

0.49 

0.82 

0.35 

GEM10 

C 

1.43 

1.75 

0.33 


A 

2.40 

2.84 

3.40 


HO 



TEST 5 DESCRIPTION: Two three day arcs are computed which overlap for one day. 
The orbits are determined from laser data. These orbits are then differenced over their 
common day and from these differences an RMS is computed for each of the ballistic 

components. 




ORBITAL COMPARISON 

TWO 3° 

ARCS WITH ONE DAY OVERLAP. 

DIFFERENCED FOR ONE DAY 


5/15 -18 vs. 18 

-21 

GRAV. 

MODEL 


3^.3° 


R 

3.07 

GEM7 

C 

7.61 


A 

41.45 


R 

0.23 

GEM9 

C 

3.19 


A 

7.14 


R 

0.32 

GEM10 

C 

3.14 


A 

7.05 


6/10 — 13 vs. 13 

- 16 


R 

4.29 

GEM7 

C 

3.93 


A 

40.73 


R 

0.60 

GEM9 

C 

1.00 


A 

5.02 


R 

0.57 

GEM10 

C 

1.02 


A 

4.94 


Ill 



TEST 6 DESCRIPTION: These tests require two steps. First, a 5 day (August 2 to 7, 1975) User 
orbit is determined with a station removed from the solution. This recovered orbit h n held 
fixed and each individual pass of this station'* User data is used to estimate a correction to the 
station height. These height (Ah) corrections arc interpreted as an estimate of the radial orbit 
error in the five day arc. This is similar to the station navigation/orbital error estimate* performed 
at NWL. 


STATION NAVIGATIONS FOR 6° ARC: 8/3 TO 8/8 

1 Station is rero weighted and orbit ; s determined. 

2 Orb«« from 1 field fixed, and station height is adjust'd for £<%ch individual pass 
of data. 


STATION 

TIME OF PASS 

RMS 

MAX. El fEV. 

AH 

crtlas 

804 0959 

1.42 

57?6 

206 

GRTLAS 

804 1935 

.37 

34°6 

0.57 

GRTLAS 

805 0944 

.82 

87°5 

—1.12 

GRTLAS 

806 204a 

.17 

68?3 

C.18 

GRTLAS 

807 2031 

.14 

75?0 

-0.02 

ARESAO 

805 0813 

2.21 

60? 2 

1.43 

ARESAO 

807 2158 

1.58 

74°4 

0.61 

ARESAO 

808 0909 

1.97 

58?3 

-0.60 

ARESAO 

808 2144 

2.05 

62°5 

1.25 

OLISAO 

S04 0144 

2.05 

58?2 

1 15 

OLISAO 

805 0130 

2.67 

80°4 

2.01 

OLISAO 

806 0116 

1.12 

56?S 

-0.57 

OLISAO 

808 0227 

1,97 

43°1 

-2.44 
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TEST 7 DESCRIPTION. 



T3 
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TEST 8 DESCRIPTION: A five day laser arc is compared with a well tracked 3 revolution arc. The 
data contained in the three revolution se9ment of the five day arc is deleted from the five day arc 
recovery. The orbital trajectory differences over this common 3 revolution vme span are statistically 
evaluated in each ballistic component as an RMS difference. 


GRAVITY MODELS DATE RMS DIFFERENCES 


5^ 

3 REV 

5^ 

3 REV 

RADIAL 

CROSS 

ALONG 

GEM9 

GEM9 

6/19-241975 

6/21 15 h - 6/21 20 h 

0.79 

0.52 

2.86 

GEM10 

GEM10 

6/19 - 24 1975 

6/21 15 h - 6/21 20 h 

0.82 

0.75 

3.02 

GEM9 

GEM9 

8/2 - 7 1975 

8/3 8 h - 8/3 14 h 

0.64 

1.31 

1.37 

GEM10 

GEM10 

8/2-7 1975 

8/3 8 h - 8/3 14 h 

0.57 

1.04 

1.20 
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TEST 9 DESCRIPTION: Three laser arcs of five days length are computed from laser data. 
This test compares the RMS of fit to the laser ranges themselves. 


RMS OF FIT (METERS) 



EPOCH 1 

EPOCH 2 

EPOCH 3 

MODEL 

750622 

750704 

750729 

GEM7 

6.84 

7.62 

11.53 

GEM8 

7.47 


6.89 

SA04.3 

12.0 2 


14.17 

GRIM2 

11.50 

11.82 

12.88 

GEM9 

1.46 

1.92 

1.25 

GEM10 

1.46 

1.91 

1.25 


;is 
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